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ABSTRACT 

The lightcurve of CoRoT-2 shows substantial rotational modulation and deformations of the planet's transit profiles 
caused by starspots. We consistently model the entire lightcurve, including both rotational modulation and transits, 
stretching over approximately 30 stellar rotations and 79 transits. The spot distribution and its evolution on the 
noneclipsed and eclipsed surface sections are presented and analyzed, making use of the high resolution achievable 
under the transit path. 

We measure the average surface brightness on the eclipsed section to be (5 ± 1) % lower than on the noneclipsed 
section. Adopting a solar spot contrast, the spot coverage on the entire surface reaches up to 19 % and a maximum 
of almost 40 % on the eclipsed section. Features under the transit path, i.e. close to the equator, rotate with a period 
close to 4.55 days. Significantly higher rotation periods are found for features on the noneclipsed section indicating 
a differential rotation of AQ > 0.1. Spotted and unspotted regions in both surface sections concentrate on preferred 
longitudes separated by roughly 180°. 
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1. Introduction 



The space-based CoRoT mission (e.g. Auvergne et al 



20091, launched in late 2006, provides stellar photome 



try of unprecedented quality. One of CoRoT's primary 
tasks is the search for extra-solar planets using the tran- 
sit method. So far several systems with eclipsing exoplan- 
ets were found, one of them CoRoT-2 harboring a giant 
close-in 'hot-jupiter'. 

The planet CoRoT-2b was discovered by |Alonso et al.| 



(2008), who determined the system parameters from tran- 
sits and follow-up radial velocity (RV) measurements. 



Bouchy et al. 



( 2008 ) observed the Rossiter-McLaughlin ef- 



fect using additional RV measurements and determined the 
projected angle A = (7.2 ± 4. 5)° between the st ellar spin 
and the planetary orbital axis. Lanza et al. (2009) used the 
strong rotational modulation of the lightcurve to study the 
spot distribution and its evolution; they detect a stellar ro- 
tation period of P rot — (4.522 ± 0.024) d and two active 
longitudes on opposite hemispheres . The seconda r y ec lipse 



of the planet was first detected by Alonso et al. (2009) in 



white light. Later the analysis was refined and extended by 



Snellen et al.| fl2009|) , detecting sign ificant thermal emission 
of the planet; Gilion et al. (2009) repeated this work us- 
ing additional infrared data. They also find an offset of the 
secondary transit timing indicating a noncircular orbit. 

The determination of the planetary parameters by 



Alonso et al. (2008) did not account for effects of stellar 
activity, which is in general not negligible for active stars. 
Czesla et al. (2009) re-analyze the transits and derive new 
parameters for the planet radius R p /R s and the orbital in- 
clination i considering the deformation of transit profiles 
due to spot occultation. This is especially important for 



attempts to reconstruct active surface regions from tran- 
sit pr ofiles. An analysis of a single transit by | Wolter et al.| 
( 2009 ) shows the potential of eclipse-mapping and yielded 



constraints on the properties of the detected starspot. 
Similar approaches to analyze signatur es of starspots in 



(2007) 



transit profiles we re also carried out b y| Pont et al^ 
(HD 189733) and |Rabus eTaTJ pO^ (TrES-1) primarily 
using data obtained with the Hu bble Space Te l escope . In 



a more comprehensive approach, Huber et al. (2009]) re 



constructed an interval of the CoRoT-2 lightcurve over two 
stellar rotations including the transits. 



This work is based on the paper of Huber et al. (2009). 
We now analyze the entire data set of CoRoT-2, modeling 
both the rotational modulation of the global lightcurve and 
transits simultaneously. In this way we derive stellar surface 
maps for both the noneclipsed and the eclipsed section of 
the star. 



2. Observations and data reduction 

The data were obtained in the first long run of the CoRoT 
satellite (May 16 to Oct. 15, 2007). The planetary system 
CoRoT-2 (Star CorotID 0101206560) consists of an active 
solar-like G7V star and a large planetary companion on a 
close orbit. Due to its edge-on view and an orbit period of 
only about 1.7 days, this lightcurve contains about 80 tran- 
sits, roughly 3 during each stellar rotation. 

The large surface inhomogeneities of this very active 
star are clearly visible in the rotational modulation of 
the CoRoT lightcurve taken over approximately 140 days. 
Although there is surface evolution even on timescales of 
one stellar rotation, these changes are small compared to 
the modulation amplitude. Significant deformations of the 
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Table 1. Stellar/planetary parameters of CoRoT-2a/b. 



Star 



Value ± Error 



Ref 



P a (4.522 ± 0.024) d L09 

Spectral type G7V B08 



Planet 



Value ± Error 



Ref. 



P p (1.7429964 ± 0.0000017) d A08 

T c [BJD] (2 454 237.53362 ± 0.00014) d A08 

i (87.7 ± 0.2)° C09 

R p /Rs (0.172 ± 0.001) C09 

a/R s (6.70 ± 0.03) A08 

u a ,u b (0.41 ±0.03), (0.06 ±0.03) A08 

P a - stellar rotation perio d 



taken from Lanza et al. 



[A08] , |Bouchy et al.| j2M§| [mEf, 



rbital 



(2009) [L091, Alonso et al 



perio 



Czesla et al] (2009) 
time of first transit, i 



(2008) 



central time ot first transit, i - or- 
bital inclination, R P ,R S - planetary and stellar radii, a - semi 
major axis of planetary orbit, u a ,u b - linear and quadratic limb 
darkening coefficients. 



transit profiles due to spots are also visible throughout the 
whole time series. 

The stellar and planetary parameters used throughout 
this analysis are given in Tab. [I] The extensive raw data 
analy s is and red uction follows the descriptions in |Czesla| 
et al.| ( |2009| ) and |Huber et al.| ( |2009[ ). 

In this paper we only analyze the alarm mode data 
of CoRoT-2, when the satellite switched from a sampling 
rate of 1/512 s _1 to 1/32 s _1 (which started after 3 tran- 
sits), and use the combination of the three color channels 
('white' light). Our analysis starts after one stellar rotation 
at JD = T c + P s = 2454242.05562 (see Tab. [I]) or a stel- 
lar rotation phase of <p s = 1.0, respectively. We use units 
of stellar rotation phase in our analysis, the last point of 
our data interval corresponds to c/) s — 31.24. This leaves us 
with an observation span of more than 30 stellar rotations 
containing 79 transits. In the following the first transit in- 
side our data interval is labeled 0, the last 78. 



3. Analysis 

3.1. The Model 

The projected axes of the planetary orbit a nd the stellar ro- 
tation are co-aligned (Bouchy et al. 2008[ ), which strongly 
suggests a 3-dimcnsional alignment. Hence, the rotation 
axis of CoRoT-2a is inclined by approximately 88°. While 
this impedes the reconstruction of latitudinal information 
of surface features, the existence of a planet crossing the 
stellar disk allows to access latitudinal information on spots 
beneath its path. During a planetary passage the surface 
brightness distribution is mapped onto the lightcurve as de- 
formations of the transit profiles. As a consequence of the 
co-aligned orientation of the planetary orbit and the stel- 
lar spin, the surface band scanned by the planetary disk 
remains the same: the planet constantly crosses the lati- 
tudinal band between 6° and 26°. Accordingly, the stellar 
surface can be subdivided into two se ctions: the eclipsed 
section and the noneclipsed section (cf., Huber et al.|2 009). 

Our surface model subdivides the two individual sec- 
tions into a number of 'strips'; N e is the number of strips 



in the eclipsed and N n in the noneclipsed section, respec- 
tively. Each strip represents a longitudinal interval inside 
the latitudinal boundaries of the corresponding section. The 
layout of our surface model is shown in Fig. 1 of |Huber et al.| 

pol. 



3.1.1. Model resolution and error estimation 

The problem of lightcurve inversion is well known to be 
ill-posed, so that the parameter space is usually further 
constraint by a regularization. One such regularization is 
the maximum entr opy approach, applied, for example, by 
Lanza et al. ( 2009 ) in their analysis of CoRoT-2a. 



We use a Nelder-Mead (N M) Simplex algorithm for min- 
imization (Press et al. 1992[ ). Our model does not require 
any regularization because o f its relatively smal l number of 
parameters. As discussed by Huber et al. (2009), we choose 
a number of strips balancing the improvement in \ 2 an d 
the deterioration in uniqueness. For higher strip numbers, 
adjacent strips increasingly influence each other because 
brightness can be redistributed without significant loss of 
fit quality. In our error analysis we assume that there is a 
unique best-fit solution to our problem and that the NM 
algorithm approaches it to within the limits of its ability 
to converge. Starting from the yet unknown best-fit solu- 
tion, brightness can be redistributed among the strips at 
the expense of fit quality. A set of Nelder-Mead fit runs 
will, therefore, provide a sequence of solutions with differ- 
ent realizations of this brightness redistribution. 

We calculate 50 reconstructions with randomized start- 
ing points and adopt the average of all reconstructions as 
our most appropriate model. As an estimate for the error 
of the strip brightness, we use the standard deviation of the 
parameter values, obtained from the set of reconstructions. 

We do point out that the averaged solution appears 
smoother than most individual reconstructions in the sense 
that the difference between adjacent strip brightnesses is 
smaller. This can be understood in the picture of brightness 
redistribution, because the brightening of one strip may 
preferably be compensated by darkening an adjacent one, 
which increases their contrast. This effect is largely can- 
celed out by the averaging, which makes it appear much 
like a regularization of the solutions. 

It may be criticized that the brightness error we use 
is largely determined by the ability of the NM algorithm 
to converge to a unique solution. We emphasize, however, 
that the x 2 range covered by the 50 reconstructions exceeds 
that required for 'classical' error analysis, and the estimate 
will, therefore, remain rather conservative from that point 
of view. 

Our test runs indicate that for our purpose the most 
appropriate number of strips to choose for the noneclipsed 
section is N n — 12, which will be used in our analysis. 
Larger numbers of N n appear to already oversample the 
surface significantly. The strip number for the eclipsed sec- 
tion is chosen to be N e = 24 which approximately reflects 
the size of feature s that are resolvable inside of transits 
( [Huber et al-1[2009] ) . 



3.2. Normalization 

The observations are normalized with respect to the ab ini- 
tio unknown spot-free photospheric flux of the star, which is 
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defined as maximum brightness b p hot = 1- Unfortunately, 
it is not trivial to obtain this photospheric flux level be- 
cause spots are likely to be located on the visible disk at 
all times. 

A possible solution for this problem is to adopt the max- 
imum observed flux as photospheric. However, this presum- 
ably introduces an error because the brightest part of the 
lightcurve shows only the flux level of the stellar disk dur- 
ing the minimum observed spot coverage, which needs not 
be zero. Lanza et al. ( 2009 ) determine an average minimum 
of flux deficit of approximately 2.5 %. Tests with slightly 
varying maximum brightness values in our reconstructions 
show no qualitative difference in the brightness distribu- 
tion except for a change in the average total spot cover- 
age; however, they show a significant decrease in x 2 of the 
lightcurve reconstructions for a maximum brightness 1 % 
to 2 % larger than the highest observed flux. As a result 
we choose a photospheric flux level of 2 % higher than the 
maximum observed flux. The entire lightcurve is normal- 
ized with respect to this value. 

3.3. Lightcurve modeling 

For analysis, the lightcurve is split into equally sized inter- 
vals, each covering one stellar rotation and three transits. 
Due to surface evolution detectable on timescales smaller 
than one stellar rotation, we choose to define a new inter- 
val after each transit. Thus, the resulting intervals overlap. 
Interval contains the transits number to 2, interval 1 
contains the transits number 1 to 3, and so on. This way 
we end up with 76 lightcurve intervals which are individu- 
ally reconstructed by our modeling algorithm. 

Using fit intervals smaller than one stellar rotation could 
further reduce the influence of surface evolution. However, 
we need at least three transits in each interval to sufficiently 
cover the eclipsed section. Therefore, we always use a com- 
plete rotation for each reconstruction interval. 

We u se th e fitting method presented in Huber 
(2009). Each fit interval is rebinned to 
= 3008 seconds for the global lightcurve and 
128 seconds inside of transits. Transit points are 



Normalized flux 



Residuals [a] 



, P P f'P'Pf'PPP , 'P'P 

~boboP^^o*o^^o*o^c^o*o 



et al. 
94x32 
to 



weighted with a factor of 10 higher than global points 
to give them approximately the same weight in the 
minimization process. We assume the planet to be a dark 
sphere without any emission; this seems to be a good 
approximation considering a secondary trans it depth of 
about (0.006 ± 0.002) % ( |Alonso et atl]2009 |. All other 
necessary parameters can be found in Table [1] 

We introduce a penalty function to suppress recon- 
structed brightnesses above the photospheric value of one. 
Without this boundary the brightness of individual strips 
exceeds this limit in some reconstructions. Strips with val- 
ues above unity must be interpreted as regions with a 
brightness greater than the (defined) photosphere; how- 
ever, in this approach we want to consider only cool sur- 
face fe atures, which was f ound to be a good approxima- 
tion by Lanza et al. (2009). Several tests showed that this 
penalty function does not significantly alter the outcome of 
our reconstructions; it primarily prevents the minimization 
process from getting stuck in (local) minima outside the 
relevant parameter range. 

To analyze the transits accurately, we require an 'undis- 
turbed transit profile', which means a profile corrected for 
the effects of stellar activity. This 'standard transit profile' 




Fig. 1. Observations and best fit model for CoRoT-2. The 
transit numbers are plotted below the lightcurve; for a de- 
tailed presentation of the transits see Fig. [9j See Sec. |3.4| 
for discussion. 



4 



K. F. Huber et al.: Planetary eclipse mapping of CoRoT-2a 



0.3 



0.2 



0.1 - 



- Z 2 r ed = 2-18 




= +0.024 
o,= 1.161 " 






H g = -0.059 
\ ° 8= L742 " 









-1 1 

Residuals [al 



Fig. 2. Distribution of residuals a between the observa- 
tions and the best fit model (see Fig.[l]). The overall distri- 
bution is the sum of two Gaussians, one coming from inside 
the transits (mean fi t and width at) and the other coming 
from outside of them (global, mean p, g and width a g ). For 
details see Sect. 13.41 



was determined by Czesla et al. ( 2009 1 , where a planet size 
of Rp/R s = 0.1 72 was found. With the planetary parame- 



ters derived by |Alonso et al. (2008) no satisfactory fits to 
the transits and global lightcurve can be produced. 



3.4. Results of the modeling 

We present our lightcurve reconstruction in Fig.[T] It shows 
the observed CoRoT-2 lightcurve (red triangles), including 
the transits, and our reconstruction (blue solid line). The 
shown reconstructed lightcurve is a combination of all mod- 
els for the 76 fit intervals; their overlaps were combined us- 
ing a Gaussian weighting. Each transit is labeled with its 
number on the lower edge of the graph; a more detailed 
picture of their fits is given in Fig. |9j where we present all 
transits in a stacked plot. 

Below the lightcurve the residuals a = (0 — C)/ao 
(O - observed lightcurve, C - reconstructed lightcurve, 
(To - error of observation) are given. The mean values for 
(To are 1.8 • 10~ 4 outside and 6.9 • 10 -4 inside of transits. 
Maximum values of the residuals are approximately 
10 a, which corresponds to a deviation of about 0.2 % 
between the observed and reconstructed lightcurves. On 
average, the absolute values of residuals are 1.12 a or 
O — C = 4.4 • 10~ 4 , respectively. 

Figure [2] shows the distribution of residuals from Fig. [T] 
It is approximately Gaussian, although the underlying dis- 
tribution is twofold: one component is a Gaussian distribu- 
tion of residuals from inside the transits with cr t = 1.161, 
and another broader component is coming from the residu- 
als of the global lightcurve, which is significantly wider with 
G g = 1.742. The addition of both Gaussians reproduces the 
overall distribution of residuals accurately. We calculate 
Xr C d = 2.18 from all residuals of the entire lightcurve. 

The interpretation of x? c d as an actual goodness-of-fit 
indicator is not straightforward in this case. An increase of 
the strip number should lead towards Xrcd values of 1. Test 
calculations show that our fit quality cannot be substan- 
tially improved beyond a certain level when the number of 
strips is increased. This level of x^ ed primarily reflects the 
evolution of the lightcurve within one stellar rotation, which 



cannot be improved within our static model, and which 
is especially visible in the global lightcurve (a g ^> at). 
Unfortunately, it is very hard to quantify this effect. 

3.5. Construction of brightness maps 

In Fig. [3] we present maps of the temporal evolution of the 
brightness distributions. The left panel gives the brightness 
map for the noneclipsed, the right panel for the eclipsed 
surface section. The rows show the reconstructed brightness 
distributions for all 76 lightcurve intervals; each interval is 
labeled by the number of the first transit it contains. Due 
to the different resolutions in the two surface sections, we 
linearly interpolate the brightness values in each individ- 
ual row; there is no interpolation applied between different 
rows. To generate the combined brightness map shown in 
Fig. [8] (left panel), the maps of the two separated sections 
are weighted, corresponding to their disk fraction of 0.79 
for the noneclipsed and 0.21 for the eclipsed section, and 
added. The errors are combined the same way. Adjacent 
rows in our maps are not independent because the fit in- 
terval is only shifted by one transit (« 1 /3 stellar rotation) 
when moving from one row to the next. 

The errors are displayed directly below each map. For 
the brightness values of the global lightcurve fit (Fig. J3j 
left panel), the mean error is about 1 %, for the transits 
(Fig. [3j right panel) it is approximately 3 %. The errors in 
the latter map are larger (on average) and more inhomoge- 
neously distributed because they also reflect the coverage of 
the eclipsed section by the transits; areas only marginally 
visible in transits cannot be reconstructed with high accu- 
racy. 

The stellar longitude scale of our maps runs backwards 
from 360° to 0°. This is due to our retrograde definition of 
the stellar longitude I compared to stellar phases <fi s ; their 
relation is I — (4>n ~ 4>s) ' 360° (with an integer stellar 
rotation number 4>n = [0s] + 1)- 

The identification of significant structures in these maps 
deserves some attention. The brightness information is color 
coded using black color for the darkest structures and white 
for photospheric brightness. Considering the approximate 
mean error of each map, we indicate areas with a brightness 
significantly below unity with yellow color. Hence, not only 
black areas of these maps are spots but yellow structures 
represent a significant decrease in brightness as well. 



4. Discussion 

In this section we discuss and interpret the spot distribu- 
tions of our brightness maps and their evolution. This in- 
volves quite a few different aspects, which are often difficult 
to disentangle, and on which we focus on individually in the 
following subsections. 

4.1. Identifying physical processes and detection limits 

Our brightness maps allow us to witness the evolution of the 
stellar surface and, to some degree, to discriminate between 
individual processes causing changes of the spot distribu- 
tion: Emergence/dissociation, differential rotation, or mi- 
gration of surface features leave potentially distinguishable 
signatures. Unfortunately, the high diversity of CoRoT-2a's 
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Fig. 3. Left panel: Brightness map reconstructed from the global lightcurve (noneclipsed surface). Right panel: Brightness 
map for the reconstructions of the transits (eclipsed surface). The combination of both maps is presented in Fig.[8j Each 
row presents the reconstructed brightness distribution of one fit interval; the transit number indicates the number of this 
interval's first transit. Each step in transit number equals a temporal step of 1.74 days. The brightness is color-coded, the 
maximum photospheric brightness is unity. The error maps at the bottom of each panel show the estimated reconstruction 
error. See Sects. [3~3l and 13.51 for details. 



surface and a limited resolution complicate the interpreta- 
tion of these signatures. 

On the eclipsed section the position of a feature is fairly 
well known and it is likely to be physically coherent. On 
the noneclipsed section it is not clear whether a feature is 
actually a single connected active region or a superposition 
of two (or even several) different regions. 

Features on the eclipsed section provide a valuable 
reference point for the detection of differential rotation. 
Systematic longitudinal movements of spots in the transit 
map do not indicate differential rotation but rather a differ- 
ence to the input rotation period. However, a comparison 
of transit and global maps could reveal different rotation 
periods on different surface sections. 

Systematic longitudinal movements of structures in the 
global map do not necessarily indicate differential rotation 
either. A simultaneous decay and growth of two distinct 
active regions at different longitudes might leave a signa- 
ture similar to a single differentially rotating active region. 
Therefore, the processes of differential rotation and spot 
evolution are hard to discern. 

A very interesting scenario is the possibility to find 
spot migration due to the planet. If a spot moves from the 
noneclipsed to the eclipsed section, its signature also moves 
from the global to the transit brightness map. If present, 
such signatures are detectable in high quality brightness 
maps. 

In the following discussion we attempt to attribute sig- 
natures in the brightness maps to the specific processes 
discussed above. Depending on the number of longitudinal 
strips used, map structures have an estimated error in lon- 
gitude of about half the strip width, which is about ±7.5° 
for the transit and ±15° for the global map, respectively. 

4.2. Spot coverage 

The brightness maps of Fig. [3] show clear evidence of (a) 
a coverage of a large fraction of the surface with dark 



features and (b) a substantial evolution of this spot dis- 
tribution within the roughly 140 days of observations. 
Considering the pronounced rotational modulation of the 
lightcurve, and the persistently changing shape between in- 
dividual rotations, this result is hardly surprising. The low- 
est flux of the rotationally-modulated lightcurve is about 
0.92 (m odulo the uncertainty in lightcurve normalization, 
see Sec. 3.2); if the starspots were absolutely dark, they 
would still cover 8 % of the disk. The largest peak-to-peak 
variation during one rotation spans from 0.98 to 0.92 indi- 
cating roughly 6 % more spots on the darker hemisphere. 

Using a spot contrast of c s = 0.7, which is about the 



average bolometric contrast of sunspots ( Beck & Chapman 



1993 



Chapman et al. 1994), we determine a maximum 



spot coverage of 37 % on the eclipsed and 16 % on the 
noneclipsed section. For the entire stellar surface a maxi- 
mum and minimum spot coverage of 19 % and 16 % are 
derived. 

Details on brightness values and spot coverage can be 
found in Fig. [4] The average brightness B is the mean 
brightness of each reconstruction interval. The top panel 
gives the ratio between the mean brightness of the eclipsed 
and noneclipsed sections, the second and third panel (from 
top) the eclipsed and noneclipsed brightnesses separately, 
and the bottom panel shows the variation of the to- 
tal brightness B Tota i = (1 — A) ■ B ecL + A ■ B non ^ ecL with 
A = 0.79. The associated spot coverage fraction is calcu- 
lated using (1 — B)/(l — c s ). 

Instead of converting brightness into spot coverage, we 
can also reverse the process. Assuming the darkest element 
of the transit map is entirely covered by one spot, the re- 
constructed brightness represents the spot contrast. One 
strip on the eclipsed section has a size of about 1 % of 
the entire surface. We obtain a value of 0.76 for the dark- 
est surface element of the transit map, which is not far 
from the solar spot contrast of 0.7. The minimum bright- 
ness of 0.76 can be translated into a temperature contrast 
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Fig. 4. Average brightness values of the eclipsed {B ec i) 
and noneclipsed {B non _ ec i.) sections for each reconstruc- 
tion interval. The y-axes on the right show the correspond- 
ing spot coverage for a spot contrast of 0.7. Top panel: Ratio 
Becl./B non -ecl. of the mean brightnesses. Second panel from 
top: Mean brightness and spot coverage for the eclipsed sec- 
tion. Third panel from top: Mean brightness and spot cover- 
age for the noneclipsed section. Bottom panel: Total bright- 
ness B rotai ', the black curve shows the results of Lanza et al. 
(2009} shifted by +9 % (see Sect. [421. 




4.5 4.75 

Period [days] 



Fig. 5. Section of the Lomb-Scargle periodogram of the 
CoRoT-2 lightcurve containing the highest peak. The three 
peaks are labeled Pi = 4.53, P 2 = 4.95, and P 3 = 4.72 days 
according to decreasing power. The adop ted rotation period 
of our reconstructions, P s — 4.522 days ( Lanza et al.|2 009), 
is close to Pi . 



of 

is at T, 



400° K betwe en the spot and the photosphere, which 
5625° K pouchy et aL|[2008| ). 



4.3. Rotation period 



For CoRoT-2a Lanza et al. (2009) determine a rota- 



tion period of (4.52 ± 0.14) d by means of the Lomb 
Scargle periodogram, which the authors later refine to be 
P s = (4.522 ± 0.024) days by minimizing the longitudinal 
migration of their active longitudes. We adopted P s for 
our reconstructions. It is almost identical to the largest 
peak in the Lomb-Scargle periodogram (Fig. p} centered 
at Pi = 4.53 d. Pi is accompanied by two other distin- 
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Fig. 6. The global surface map (left panel of Fig. i| 
with pronounced dark and bright features marked and la- 
beled Dl - D3 and Bl - B5, respectively. The (green) dash- 
dotted line in dicat es the tentative movement of bright re- 
gions (see Sec. 4.4). 



guishable peaks at P 2 = 4.95 d and P3 = 4.72 d, which are 
related to structures in our surface maps as discussed in 
Sect. l44l 

Since the planet crosses the stellar disk in a latitudinal 
band between 6° and 26°, spots on the eclipsed section must 
be located close to the equator. The approximate vertical 
alignment of the darkest features in the transit map indi- 
cates that low-latitude features rotate close to the adopted 
rotation period of 4.522 days used in our reconstructions. In 
contrast, features on the noneclipsed section do show longi- 
tudinal migration, which may be attributed to differential 
rotation of spot groups at latitudes > 30°. This finding sug- 
gests that active regions close to the equator dominate the 
modulation of the lightcurve. 

A closer inspection of the transit map shows a small 
but constant longitudinal shift of about 60° over the en- 
tire 30 rotations. It is not only visible in the dark struc- 
tures located at ~ 200° and ~ 60° longitude, but for the 
bright region at ~ 300° as well. As a consequence, these 
low-latitude features do not rotate exactly with a rotation 
period of 4.522 days but about 2° per rotation more slowly, 
which tra nslates into a rotat ion period of approximately 
4.55 days. Lanza et al. (20091 state a retrograde migration 
of their second active longitude corresponding to this rota- 
tion period. 

We recalculated the surface reconstructions adopting a 
rotation period of 4.55 days. As expected the longitudinal 
shift previously detected in the transit map disappeared, 
but some features now seem to migrate in the other di- 
rection. Therefore, the exact rotation period of features 
on the eclipsed section is probably slightly smaller than 
4.55 days. The global map changes accordingly when using 
a rotation period of 4.55 days; the brightness distribution 
reconstructed at larger transit numbers is shifted towards 
larger longitudes, but remains qualitatively the same, so 
that these maps are not presented separately. 

4.4. Longitudinal movement / differential rotation 

Some attributes of the global map suggest a longitudinal 
movement of surface structures; a labeling of particularly 
interesting regions is given in Fig. [6] Especially considering 
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the inactive regions, the bright structure Bl at the bottom 
left appears to move from ~ 300° to ~ 100° (B2) within 
about 4 stellar rotations. It continues moving in the di- 
rection of decreasing longitude reaching roughly 240° near 
transit number 50 (B3). Regions B4 and B5 do not seem to 
fit well into the line drawn by Bl, B2, and B3. However, a 
line connecting B4 and B5 roughly matches with the pre- 
viously detected rotation period between 4.8 and 5.0 days. 
The tilted shape of the dark structure Dl roughly agrees 
with this range of periods as well. D2 indicates a smaller 
rotation period consistent with P3. 

Although such surface map characteristics do not nec- 
essarily prove differential rotation, they are certainly sug- 
gestive of differential rotation. Assuming now this to be 
correct, we will elaborate on its consequences. The largest 
rotation periods we obtain from structures of the global 
map lie between about 4.8 and 5.0 days. We studied the 
longitudinal movements of these structures on surface maps 
obtained from lightcurve reconstructions applying different 
rotation periods and always end up with similar results. An 
examination of CoRoT-2's periodogram reveals a splitting 
up of the highest peak into three components: P\ = 4.53 d, 
P 2 = 4.95 d, and P 3 = 4.72 d (see Pig.|5|. It is striking that 
P2 is close to the rotation period determined from tilted 
structures in our global map. 

If these rotation periods do arise from differential rota- 
tion, we can estimate a lower limit of its strength. The 
rotation period of eclipsed spots close to the equator is 
consistent with the highest peak P\ = 4.53 days of the peri- 
odogram, the largest periods detected in the global map are 
around P 2 = 4.95 days; therefore, we adopt this period for 
the most slowly rotating active regions. This way we deter- 
mine a lower limit of AO > 0.1 rad/d or a > 0.08. This is 
consistent with values expected for stars with temp eratures 
and rotatio n periods similar to that of Co RoT-2 (Barnes 
|et al.|2005| . Using a 3-spot model approach, Frohlich et a 
(2009J derive an estimate of AO > 0.11 rad/d, which is in 
good agreement with our result. 

Although the peaks in the periodogram ht nicely in with 
our brightness maps, attributing them to three active re- 
gions with associated rotation periods may not be fully ad- 
equate. We simulated several lightcurves with differentially 
rotating spots and examined their periodograms. Although 
the main peak splits up into different components, the peri- 
odogram does not necessarily map the exact rotation peri- 
ods of the differentially rotating spots to the peak barycen- 
ters. An exhaustive analysis of the periodogram is beyond 
this work's scope, but the above approach may serve as an 
approximation, and it shows that the characteristics of the 
periodogram can be aligned with surface map attributes. 

As an alternative, or maybe extension, to the interpre- 
tation in terms of differential rotation, an evolution of the 
global activity pattern should be considered, which does 
not invoke longitudinal movement of individual surface fea- 
tures, but a redistribution of strength between active re- 
gions. The spot distribution on the global map suggests 
sudden longitudinal relocations of the most active feature. 
For the first ~ 25 transits the dominant spotted region 
Dl keeps its position at ~ 30° showing an apparent move- 
ment towards smaller longitudes at the end of this interval. 
Afterwards the dominant active region is found approxi- 
mately 180° apart from the previous position in region D2. 
Such 180°-jumps, or rather 'hemisphere-jumps', as the value 




lude [degrees] 



Fig. 7. Direct comparison of the two brightness maps of 
Fig. [3] The global map is drawn in color, structures of the 
transit map are shown with contour lines. Upper panel: Map 
of bright structures in the global map with b non - ec i, > 0.95. 
The contour line delimits brightness values of the transit 
map above 6 ec /. = 0.94. Lower panel: Map of dark struc- 
tures with b non —eol. < 0.95. Here the contour lines indicate 
a transit map brightness of b ec i, = 0.87. 



of 180° should not be taken too seriously here, are found 
for inactive (bright) regions as well. 

An appealing explanation for this apparent shift is a 
change of the relative strength of two active regions, which 
does not involve movement of any of the structures them- 
selves. These jumps are possibly a sign for s ome flip-flop 
scenario as already claimed for other stars (Jetsu et al. 
1993 Korhonen et al. 2001), where the relative strength 
between two active longitudes is changing suddenly and 
periodically. The timescale of these 'jump periods' derived 
from our global brightness map is roughly 10 stellar rota- 
tions; however, this is not seen in the transit map. 



4.5. Lifetimes of features 

Both the global and the transit map provide the possibility 
to measure lifetimes of spotted regions. For the noneclipsed 
section it is not clear whether a dark structure actually 
is a single connected active region or a superposition of 
several individual ones at roughly the same longitude but 
completely different latitudes. 

The transit map shows a high stability of the spot po- 
sitions. One group of spots is located between 200° and 
240°, which is stable up to transit number 35; then the spot 
coverage seems to decrease for about one stellar rotation. 
Afterwards the spot distribution becomes more complicated 
and spreads over a larger area. The left part of this more 
complex region, at about 210° longitude, is probably a con- 
tinuation of the preceding active region; however, at about 
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160° longitude a new, clearly separated group of spots ap- 
pears. This latter structure appears to have a lifetime of 
approximately 10 stellar rotations, which would be in good 
agreement with timescales observed on the global map. 

Another group of spots is located between 0° and 100°. 
At the beginning of the observation these spots are spread 
over this larger interval of about 100° on the eclipsed sec- 
tion; later, approximately after transit number ~ 50, they 
seem to concentrate in a smaller interval. This active re- 
gion seems to have a significant fine structure in time, 
which might indicate smaller lifetimes of individual spot 
groups; however, those structures could also be artifacts of 
the reconstruction algorithm. We prefer the interpretation 
that this active region, and the left side of the other ac- 
tive region, are persistent over all 30 rotations, although 
they probably contain smaller spots undergoing significant 
evolution on much shorter timescales. 

The left side of both maps (around 300° in longitude) 
shows bright regions which are least covered by spots during 
all observations, indicating a lifetime of half a year for this 
'inactive longitude'. For the darkest structures of the global 
map, maximum lifetimes of about 10 to 15 stellar rotations 
can be estimated. This is in agreement with the results of 
Lanza et al. (2009), who determine a lifetime of 55 days 
(=12 stellar rotations) for active regions, and identify this 
time span with the 'beat period' visible in the lightcurve. 

Although there are some structures in the transit map 
that suggest smaller lifetimes than half a year for individual 
spot groups, the lifetimes of features on the two different 
surface sections seem to be different. A direct comparison of 
the global and the transit map is shown in Fig. [7] The active 
regions on the eclipsed section remain active all the time 
despite their possibly significant fine structure. In contrast, 
active regions on the noneclipsed section evolve faster show- 
ing more pronounced changes and a longitudinal movement 
compared to the spots on the eclipsed section. The reason 
for the apparent lifetime difference is not clear. Possibly the 
darkest structures represent only a superposition of several 
spotted groups at about the same latitude. If these groups 
change their mutual longitudinal positions, or a fraction of 
the spot groups dissolves, the dark structures in our maps 
would brighten. In this case the darkest structures would 
only represent special configurations of the spot distribu- 
tion and their 'lifetimes' in our maps would not be directly 
connected to the lifetimes of individual active regions on 
the surface. 

The stable vertical alignment of features in our transit 
map cannot be caused by a systematically incorrect tran- 
sit profile used in our reconstructions. The transits do not 
always cover exactly the same part of the eclipsed section, 
which is also visible in the error distribution of the transit 
map; dark structures indicate where the coverage is best, 
bright where it is worst. Thus, an error introduced by the 
transit profile would be distributed over the entire map. 

4.6. Comparing global and transit maps: spot migration? 

Finally, there is the possibility of detecting signatures of 
spot migration - the movement of spots from the equator 
to the poles or vice- versa - in our brightness maps. Figure]?] 
presents the direct comparison of the brightest structures 
(upper panel) and the darkest structures (lower panel) of 
the global and the transit map. Especially the dark struc- 
ture D2 (see Fig. |6|) in the middle of the lower panel's map 



suggests that features on the eclipsed and noneclipsed sec- 
tions are related. For the first 30 transits there is a dark 
feature of the transit map at this longitude, then it starts 
to disappear when D2 becomes darker. This reflects a sce- 
nario where a spot group moves from the eclipsed section 
to the noneclipsed. After transit number 40, D2 starts to 
disappear while other structures appear on the eclipsed sec- 
tion. If this really represents a case of spot migration, the 
spot group either moves back onto the eclipsed section, or it 
stays outside and new spots emerge under the transit path. 

A similar observation can be made concerning the bright 
structures in the upper panel of Fig. [7| The bright struc- 
tures in the transit map between 60° and 120° alternate 
with the bright regions B2 and B4. First there is a bright 
structure on the eclipsed map below region B2, then there 
is a little bit of both between B2 and B4, and after region 
B4 a bright structure is emerging in the transit map. 

It is impossible to prove whether these signatures really 
represent spot migration; probably some of them are due 
to other processes, e.g. short-term evolution of spotted re- 
gions. Nevertheless, there is a similarity to what one would 
expect to see in brightness maps from surfaces showing spot 
migration. A behavior supporting a shift of spots from the 
eclipsed to the noneclipsed sections (and vice-versa) can be 
observed in Fig. [4] (second and third panel). It suggests a 
correlation between the mean brightnesses of the two sec- 
tions; when the average brightness of the eclipsed section 
decreases, it increases on the noneclipsed part. However, 
this correlation does not necessarily prove a steady motion 
between the two sections and might as well indicate that 
vanishing spots just reappear somewhere else. 

4. 7. Comparison to previous results 

Figure [8] (left panel) displays the combined brightness 
map derived from both the eclipsed (transit map) and 
noneclipsed sections (global map) of Fig. [3] the single maps 
are multiplied by their corresponding surface fractions (0.21 
for the eclipsed and 0.79 for the noneclipsed) and added. 



Lanza et al. ( 2009 ) present a map of the surface evolu- 
tion derived from a fit to the global lightcurve (their Fig. 4) 
not including the transits. In the right panel of Fig. [8] we 
present a comparison of their results to ours. Since we do 
not use filling factors, we translated their map into bright- 
nesses using their spot contrast of 0.665. We take the result- 
ing map (color coding) and superimpose it on our combined 
map from both the eclipsed and noneclipsed sections (con- 
tours). In the left panel of Fig.|8]the same contour lines are 
drawn to provide a better comparison. |Lanza et al.f s and 
our results show good agreement, although a perfect match 
in fine-structure is neither found nor expected. Dark and 
bright structures are located at very similar positions and 
the shapes are consistent. 

Adding up the brightness values of each reconstruction 
interval of the map in Fig. H (left panel), we can study 
the variations of the mean total_ brightness B^otai of the 
star. This is presented in Fig. 
maximum of Bxotai — 0.951 anc 
the maximum difference between the highest and lowest 
average total brightness is only about 1 %, which is much 
less than the maximum brightness differences within the 
brightness maps. This implies that the star as a whole does 
not change its overall spot coverage as dramatically as it 
redistributes it; when spots disappear, other spots show 



4] (bottom panel). With a 
a minimum value of 0.942, 
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Fig. 8. Left panel: Brightness distribution of the entire surface for all reconstructed intervals; it represents a combination 
of both the global and the transit map. The error for each bin of the brigh tness map is shown below. Right panel: 
Comparison of our brightness distribution (contours) to the reconstruction of Lanza et al. (2009), which are shown in 
color coding. See Sect. |4.7| for details. 



up. The solid line in the pan el gives the comparison to the 
results of Lanza et al. ( 2009 1 . We translated their values to 
our spot contrast of 0.7 and shifted it by a constant spot 
coverage of +9 % to match our points. Our average spot 
coverage of 17.5 % is roughly twice as high. Although, a 9 % 
shift seems to be enormous, about 80 % of it (+7 %) can 
be attributed to a different n ormalization of the lightcurve 
and, thus, photosphere. While Lanza et al. ( 2009 1 define the 
maximum flux in the lightcurve as their phot ospheric level, 
our photosphere is 2 % brighter (cf. Sect. 3.2 1, which has to 
be compensated by spots. We attribute the remaining 2 % 
to the differences in the adopted models. In particular, we 
use longitudinal strips and spots can only be distr i buted 



(2009) 



homogeneously across a strips, while Lanza et al. 
localize the spots in 200 bins on the surface 

Previously, we detected an average brightness under the 
eclipsed section (6 ± 1) % higher than on the noneclipsed 



section ( Huber et al.|2009 ). This value is redetermined from 
the reconstruction of the entire lightcurve presented in this 
paper. It decreases to (5.4±0.9) % (see top panel of FigQ. 



4.8. Brightness maps and lightcurve modulation 

It is striking that the rotational variations of the star are 
additionally modulated with a beat period about a factor 
10 to 15 larger. During the maxima of this large-scale mod- 
ulation (at about transit numbers 15, 45, and at the end 
of the lightcurve), the minima of the stellar rotation are 
deep and regular, i.e., the lightcurve has only one distinct 
minimum per rotation. During the minima of this beat pe- 
riod (at about transit numbers 30 and 60), the rotational 
modulation is flatter and more complex; the minima are 
split up in two. The beat period maxima indicate the ex- 
istence of one large active region or longitude dominating 
the stellar surface. The minima indicate that two smaller 
active regions at significantly different longitudes imprint 
their signatures onto the lightcurve leaving double-peaked 
structures during one stellar rotation. This means dark re- 
gions are redistributed on timescales of 10 to 15 stellar ro- 
tations from essentially one large feature to at least two 



smaller, longitudinally separated ones, and then back to a 
large one. 

This is also observable in our brightness maps. Figure [8] 
shows a change of the dominant surface feature from ~ 60° 
to ~ 220° at about transit number 20. Earlier the lightcurve 
was dominated by one large active longitude leading to one 
broad minimum during each rotation. During the transition 
phase, especially around transit number 25, the lightcurve 
minima become double-peaked. After transit number 30 
the activity center of the active longitude at 220° moves 
to about 180° and, thus, closer to the other active longi- 
tude at about 60°, which has not entirely disappeared and 
starts to grow stronger again. Because the active longitudes 
are closer now, the minima grow deeper; additionally, the 
bright region at 300° becomes larger and the maxima rise. 

Again there are double-peak structures between transit 
numbers 55 and 60 for the same reasons. The decrease of the 
strong maxima is primarily due to the temporary change 
of the bright longitude at 300°. Interestingly the amplitude 
and width of the minima become even larger at the end of 
the lightcurve, where the active longitudes first move closer 
together. In the end, the feature at ~ 20° becomes dominant 
spreading over almost 60° in longitude and leading to very 
low lightcurve minima. 

The main reason why the rotation period of the star is 
so nicely indicated by almost constantly separated maxima 
of the lightcurve - despite the pronounced surface evolu- 
tion - is that the on average brightest part of the surface, 
located on the inactive longitude at 300°, remains stable at 
its position. 



5. Conclusions 

We present a reconstruction of the complete CoRoT-2 
lightcurve - including transits - covering about 140 days 
or 30 stellar rotations. In contrast to previous work, both 
the transit profiles and the rotationally-modulatcd global 
lightcurve are fitted simultaneously leading to a consistent 
solution for the entire lightcurve. From the transits the 
brightness distribution on the eclipsed surface section is re- 
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Fig. 9. Observation (blue triangles) and reconstruction (red solid line) of 79 transits from the CoRoT-2 lightcurve. The 
first transit of each column is shifted to a continuum of zero, each subsequent transit is shifted by —0.025. The number 
of each transit is annotated inside the plot. See Fig. [TJfor the entire lightcurve. 



covered, which is the part of the surface between 6° and 26° 
latitude constantly eclipsed by the planet. The noneclipsed 
section is reconstructed from the rotational modulation. 

The evolution of the spot distributions on both sur- 
face sections is presented in two maps showing the surface 
brightness distribution as a function of time. The compos- 
ite map, which shows the evolution of the entire surface, 
is juxtapositioned with previously published results. The 
results are found to be in agreement taking into account 
different modeling approaches and assumptions. 



The transit map shows two preferred longitudes densely 
covered with spots and separated by approximately 180°. 
Both active regions persist for the entire observing time 
of ~ 140 days, although they undergo significant evolution 
in size, structure, and brightness. In the transit map they 
show a constant retrograde movement indicating that the 
adopted stellar rotation period of P s — 4.522 days does not 
exactly describe their rotation. We determine that these 
low-latitude features rotate at a period of approximately 
4.55 days, which is also true for the long-lived inactive lon- 
gitude at ~ 300°. 
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The global map is more complex than the transit map 
presumably because the structures it describes represent 
a superposition of features at similar longitudes but dif- 
ferent latitudes. Usually it shows only one dominant dark 
feature at a time, which changes position after approxi- 
mately 10 to 15 stellar rotations. Again these dominant 
features are separated by about 180° in longitude. A per- 
sistent inactive longitude exists at about 300° similar to 
the one in the transit map and at about the same position. 
The global map indicates that there are features located 
on the noneclipsed section with significantly larger rota- 
tion periods than 4.522 days. This suggests the presence of 
differential rotation with spots moving more slowly at high- 
latitudes than at low-latitudes. We estimate a differential 
rotation of A£l > 0.1 or a > 0.08, respectively. 

Assuming a spot contrast of 0.7, the spot coverage of the 
eclipsed section reaches a maximum of 37 %, which is more 
than twice as large as the maximum on the noneclipsed sec- 
tion. On average the eclipsed section is (5±1) % darker than 
its noneclipsed counterpart. Sunspots are located within 
±30° around the solar equator. Similarly, our results indi- 
cate that spot groups on CoRoT-2 are also concentrated in 
a low-latitude 'active belt'. 
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